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ABSTRACT 

We explore the relationships between the Polycyclic Aromatic Hydrocarbon (PAH) feature strengths, 
mid-infrared continuum luminosities, far-infrared spectral slopes, optical spectroscopic classifications, 
and silicate optical depths within a sample of 107 ULIRGs observed with the Infrared Spectrograph 
on the Spitzer Space Telescope. The detected 6.2 fim PAH equivalent widths (EQWs) in the sample 
span more than two orders of magnitude (~0. 006-0. 8 /im), and ULIRGs with Hll-like optical spectra 
or steep far-infrared spectral slopes (>S'25/S'6o < 0.2) typically have 6.2 /im PAH EQWs that are 
half that of lower-luminosity starbursts. A significant fraction (~40-60%) of Hll-like, LINER-like, 
and cold ULIRGs have very weak PAH EQWs. Many of these ULIRGs also have large (rg.7 > 2.3) 
silicate optical depths. The far-infrared spectral slope is strongly correlated with PAH EQW, but not 
with silicate optical depth. In addition, the PAH EQW decreases with increasing rest-frame 24 fim 
luminosity. We argue that this trend results primarily from dilution of the PAH EQW by continuum 
emission from dust heated by a compact central source, probably an AGN. High luminosity, high- 
redshift sources studied with Spitzer appear to have a much larger range in PAH EQW than seen 
in local ULIRGs, which is consistent with extremely luminous starburst systems being absent at low 
redshift, but present at early epochs. 

Subject headings: infrared:galaxies - galaxies: active - galaxies:starburst 



1. INTRODUCTION 

Ultraluminous Infrared Galaxies (ULIRGs) have bolo- 
metric luminosities comparable to quasars (Lboi~10^^ 
Lq), but emit nearly all of this energy at mid- 
and far-infrared wavelengths. Although ULIRGs are 
rare in the local Universe and cor nprise only a few 
percent of in frared-bright galaxies (jSoifer et al.l 1 19871 : 
ISanders et al.ll2b03> ). they account for a rapidly increas- 
ing fraction of all star- formation activity at high-redshift, 
and they may domin ate the f a r-infr a red background 
at z > 2 iFr ances chini et all [2OOII : iLe Floc'h et all 
120051: iPerez-Gonzalez et al.l 1200511 It has been sug- 
gested that they play a role in t he formation of both 
quasars and elliptical galaxies (iKormendv fc Sanders! 
fl99l ISanders fc Mirabe]lil996HScott et al.ll2002D . 

To understand the roles that ULIRGs play in the 
formation of massive galaxies, the global star for- 
mation history of the universe, and the creation of 
the cosmic infrared background, it is essential to de- 
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termine the extent to which Active Galactic Nuclei 
(AGN) and star formation contribute to their enor- 
mous luminosities. Seve r a.1 studies have used op- 
tical (dc Grii p et al. 

19951 : 
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infrared (iGoldader et al.l 119951; IVeilleux et al.l 



1999a^ and near- 
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Davies et al. |2003[ iDannerbauer et al.l l2005[ r spec 
troscopy to determine the main source of ionization of 
the line-emitting gas in local ULIRGs. These studies 
show that while the majority of ULIRGs have optical 
and near-infrared spectra similar to those of starburst 
galaxies, the fraction of ULIRGs displaying the spectro- 
scopic signatures of an AGN, either broad permitted lines 
or high-ionization narrow-lines, increases with infrared 
luminosity. In addition, evidence for AGN activity is 
more commonly found among ULIRGs with flatter far- 
infrared spectral slopes (those classified as "warm" with 
'5'25/'S'6o > 0.2 as measured with IRAS), than among 
ULIRGs with steeper far-infrared spectral slopes (those 
classified as "cold" with S25/Seo < 0.2). 

A principle uncertainty of these optical and near- 
infrared studies is the extent to which the physical con- 
ditions of the visible gas are linked to those in the deeply 
embe dded nucleus, whe re the bulk of the energy is gener- 
ated (jSoifer et al.ll2000f ). Emission at longer wavelengths 
can penetrate more heavily obscured regions close to 
the nuclear power source. Therefore, diagnostics derived 
from mid-infrared spectroscopy are potentially a better 
probe of the dominant energy source in the nuclear re- 
gions. The mid-infrared spectra of ULIRGs are com- 
prised of continuum emission from heated dust grains; 
broad emission features associated with Polycyclic Aro- 
matic Hydrocarbons (PAHs), the most prominent of 
which appear at 6.2, 7.7, 8.6, 11.3, and 12.7 fim; ah- 
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sorption by amorphous silicates centered at 9.7 and 18 
^m; and atomic fine-structure lines of Ne, O, Si, and 
S covering a large range in ionization potential. Di- 
agnostic diagrams using combinations of fine-structure 
line ratios, mid-infrared spectral slope, and PAH fea- 
ture strengths have been used to classify bright ULIRGs, 
based on the expectations that ULIRGs with central 
AGN produce high ionization lines, have flatter spectral 
slopes, and display lower PAH equivalent widths than 
those without fe.g. iGenzel et al.|[l99l iLutz et al.lll998l: 
Rigqpoulou et all Il999l: iLaurent erani200ClHTran et all 



20nit iSturm ei^i 120021 : [Armus et al.N2n04 120061 . 12007t 



Farrah et al.ll2007[ ). For sources at high-redshift (z > 1), 



however, it is generally not possible to place useful lim- 
its on the high- ionization fine-structure lines. The PAH 
emission, together with the infrared spectral slope, are of- 
ten the only mid-infrared diagnostics available for study- 
ing high-redshift ULIRGs. 

The first attempts to classify the energy generation 
mechanism in ULIRGs using mid-infrared PAH fea- 
tures m ade use of data from t h e Infrared Space Obser- 
vatory (iLutz et all Il998l 1X9991 : iRigopoulou et all Il999l : 
iTran et al.ll2'00li r For these studies, ULIRGs were clas- 
sified based on the strengths of their 7.7 /im PAH line-to- 
continuum ratios (L/C, defined as the ratio of the peak 
height of the 7.7 /im PAH feature to the level of the 
underlying 7.7 /xm continuum). The PAH-based classifi- 
cations showed that most (~80%) ULIRGs have strong 
PAH emission indicative of star formation as the dom- 
inant power source. However, the highest luminosity 
sources tend to have weak (or undetected) PAH emission, 
interpreted as a sign o f increased AGN activity. Recently, 
llmanishi et al.l (|2007| ) analyzed the mid-infrared spectra 
of 48 nearby U LIRGs observed w ith the Infrared Spec- 
trograph fIRS: iHouck et al]|2004f ) on board the Spitzer 
Space Telescope. These authors suggest that 30-50% of 
these ULIRGs, all of which have HH-like or LINER-like 
optical classifications, harbor buried AGN. 

In this paper we present the mid-infrared spectroscopic 
properties of a large sample of 107 ULIRGs that comprise 
the IRS GTO team ULIRG survey. The sample includes 
a large fraction (^40%) of warm sources, which tend 
to have higher infrared luminosities than cold ULIRGs. 
The warm ULIRGs are a critical population to study, 
since it has been suggested that they mark the tran- 
sition between staburst-do minated ULIRGs and QSOs 
(e.g. ISanders et al.l fl9"88b[ ). Because of their high lu- 
minosities, they provide a much-needed baseline over 
which to extrapolate to the extremely high luminosities 
(LiR > 10"'^^ Lq) seen at higher redshift. 

The sensitivity of the IRS offers two distinct advan- 
tages over previous studies with ISO. First, we can mea- 
sure the PAH strengths with respect to the dust con- 
tinuum (the equivalent width, or hereafter EQW) over 
nearly two orders of magnitude, allowing us to quantify 
the PAH emission in many ULIRGs that previously only 
had upper limits. This is critical for fitting trends with 
luminosity and spectral slope, as many ULIRGs have 
much weaker PAH emission features than are typically 
found in low- luminosity starburst galaxies (see H3.1^ . 
Second, the Short-Low and Long-Low IRS modules (5- 
38 /im) offer greatly increased wavelength coverage com- 
pared to ISOPHOT-S, allowing us to accurately fit the 
silicate absorption (at 9.7 and 18 /im) and multiple PAH 



features in all the ULIRGs. Here, we have chosen to 
fit the 9.7 /im silicate absorption feature since it is the 
strongest and the 6.2 and 11.3 //m PAH features be- 
cause they are both well-isolated and easily measured. 
Previously the 7.7 /im feature was often used, which is 
not ideal because absorption at both short er and longer 
wave lengths can mimic a peak at 7.7 /^m (jSpoon et al.l 
I2OO2I ) and its integrated flux can be difficult to measure 
due to blending with the adjacent 8.6 /im feature. By 
correlating the PAH emission with the silicate absorp- 
tion, it is possible to gain insight into the distribution 
of the heating source(s) and the grains. In the following 
analysis, we use Hq = 70 km Mpc^^, fim = 0.3, and 
A = 0.70. 

2. OBSERVATIONS & DATA REDUCTION 

The 107 ULIRGs presented in this paper were chose n 
primarily from the IRAS 1-Jy survey (Kim et al.lll99"8l ). 
the IRAS 2-Jy survey (jStrauss et all I1992D. and the 
FIRST /IRAS radio- far- infrared sample of lStanford et al.l 
((200'oi) . The ULIRGs span a redshift range of 0.018 < 
z < 0.93, have IRAS 60 /im flux densities between 
0.14 < Seo/Jy < 103, and have integrated 40-500 /im 
luminosities in the range 11.7 < log(LFiR/L0) < 13.13. 
The sample consists of 41 warm ULIRGs and 66 cold 
ULIRGs, based upon their rest-frame far-infrared flux 
densities at 25 and 60 /im. Man y of our highest lumi - 
nosity sources were taken from the [Stanford et all ()2000f ) 
FIRST/IRAS survey. While radio-bright AGN may pref- 
erentially populate the highest luminosities of the Stan- 
ford et al. survey, the selected ULIRGs have radio to 
far-infrared flux density ratios consistent with infrared- 
selected starburst galaxies of lower luminosity. The sam- 
ple is specifically designed to include a larger fraction 
of warm and infrared luminous sources than would be 
found in a pure flux- or volume-limited sample, for the 
express purpose of allowing a careful study of the mid- 
infrared spectral properties as functions of infrared color 
and luminosity. 

All ULIRGs were observed in Staring Mode with both 
sub-slits (orders) of each of the Short-Low (SL) and 
Long-Low (LL) modules of the IRS. Each target was ac- 
quired by performing a high accuracy IRS peak-up on the 
target itself, or by peaking up on a nearby 2MASS star 
and offsetting to the target. Each galaxy was observed at 
two nod positions within each of the IRS sub-slits (SLl, 
SL2, LLl, LL2). The resulting spectra have a spectral 
resolution of R ~ 80 over the 5-38 /im wavelength range. 

All spectra were reduced using the S14 IRS pipehne 
at the Spitzer Science Center. This reduction includes 
ramp fitting, dark sky subtraction, droop correction, lin- 
earity correction, and wavelength and flux calibration. 
The flux calibration sources were HD 1735 11 for the SL 
and LL2 modules and KsiDra for the LLl module. For a 
given module, order, and nod position, the background 
in the two-dimensional spectrum was removed by sub- 
tracting the combined two-dimensional data taken with 
the same module, but adjacent order. One-dimensional 
spectra were extracted from the background-subtracted 
two-dimensio nal spectra using th e SMART data reduc- 
tion package ([Higdon et al.l [200^ . The adopted extrac- 
tion apertures arc four pixels at the blue end of each 
order and expand linearly with wavelength. 

The 6.2 and 11.3 /im PAH features were measured by 



3 



integrating the flux above a spline-interpolated contin- 
uum. For those sources with strong 6 /im water ice ab- 
sorption, the apparent 6.2 fim EQW has been corrected 
by using an inferred 6.2 continuum defined by a spline 
interpolation between 5-26 ^m, using pivot points at 5.2, 
5.6, 7.8 (in the case of very weak PAH emission), 14, and 
26 fim. This correction results in a lower 6.2 /xm PAH 
EQW, and is correct under the assumption that the PAH 
emission is unaffected by the bulk of the ice absorption. 
A full description o f the ice-fitting procedure is given by 
ISpoon et al.l l|2007t ). who first present the measured 6.2 
fiui EQWs for many of these ULIRGs (see §3.2 and Fig- 
ure 2 for a description of the correction and Figure 1 for 
an indication of the size of the correction). The PAH 
emission at 11.3 /im is affected by the broad 9.7 /im sil- 
icate absorption feature, which can be deep in ULIRGs. 
However, lacking a strong understanding of how the ob- 
scuring dust is distributed, we have not attempted to 
correct for it. Of the 107 ULIRGs analyzed in this pa- 
per, the 6.2 /im PAH feature is detected for all but 11 
sources, and the 11.3 fim PAH feature is detected for all 
but 12 sources. 

In the following, we investigate the variation of the 
PAH features with rest-frame continuum fluxes at 5.5, 
24, 25, and 60 /im. The 5.5 /im flux is taken to be the 
average continuum betweeen 5.3 and 5.8 /im. The 24 /tm 
flux was calculated by convolving the IRS spectrum with 
the MIPS 24 /tm bandpass. In the 10 cases where the IRS 
spectrum did not extend over the full wavelength range 
of the wide IRS bandpass due to trimming of pipeline 
artifacts (e.g. fringes), it was extended by linearly in- 
terpolating in log-log space betwee n the IRS spe ctrum 
and the 60 /im IRAS flux density (iMoshid [1991 . No 
color-correction was made to place these filter fluxes on 
the MIPS absolute scale, but these corrections are small 
(<5%) based on the spectral shape of the ULIRGs within 
the 24 /tm bandpass (see the MIPS Data Handbook). 
The 25 /tm rest-frame flux was estimated by taking the 
median of the trimmed and (when needed) interpolated 
IRS spectrum between 24.7 and 25.3 /tm. The rest-frame 
60 /im luminosity was estimated by linearly interpolat- 
ing in log-log space between the observed IRAS 60 and 
100 /tm flux densities. For 12 ULIRGs, this interpola- 
tion likely underestimates the 60 /tm luminosity, either 
because the object is at a redshift high enough {z > 0.65) 
that the observed 100 /tm IRAS observation corresponds 
to a rest-frame wavelength less than 60 /tm, or because 
the 100 /tm observation resulted in only a limit, rather 
than a detection. 

Classifications based on optical spectra are available 
from the literature for 64 of the 107 ULIRGs in our sam- 
ple. Of these, 15 have Hll-like optical spectra, 22 are 
LINER-like, and 27 are Seyfert-like. 

3. RESULTS 

3.1. PAH strength and silicate absorption versus 
spectral slope and optical classification 

The median IRS spectra for warm and cold ULIRGs 
are shown in the left-hand panel of Figure [T] Cold 
ULIRGs have a median 6.2 /tm PAH EQW that is a 
factor of six higher than the median measured for warm 
ULIRGs (0.24 versus 0.04 /im; see Table [T]). The right- 
hand panel of Figure [T] shows t he average starburst spec- 
trum from iBrandl et all ()2006D plotted over the median 



cold ULIRG spectrum. Cold ULIRGs have a median 6.2 
/tm PAH EQW which is about h alf the median f o r lowe r- 
luminosity starbursts (see also llmanishi et al.l (j2007f )). 
The difference in the PAH EQWs in the median warm 
and cold ULIRG spectra reflects the strong trend be- 
tween PAH EQW and mid-infrared spectral type shown 
in Figure [21 There is a huge range in measured PAH 
EQW among the sample ULIRGs. The full range in 
PAH EQW for those ULIRGs with detected emission is 
0.006-0.864 /im for the 6.2 /xm feature, and 0.006-1.169 
/tm for the 11.3 /tm feature. The smallest upper lim- 
its are 0.005 and 0.004 /tm for the 6.2 and 11.3 /tm fea- 
tures, respectively. While the median 6.2 /tm PAH EQW 
among cold/HII ULIRGs is about half as large as for lo- 
cal starburst galaxies, the median 11.3 /im PAH EQW in 
cold/HII ULIRGs is 80% of that measured in starburst 
galaxies. The cold ULIRGs and starburst galaxies track 
each other over the limited range of 0.06 < S'25/<S'6o < 0.2 
(see Figure [2]). 

Figure[T]also demonstrates that the median cold source 
shows stronger silicate absorption at both 9.7 and 18 /tm 
than the median warm source (1.36 versus 0.65 for rg.7, 
the optical depth at 9.7 /tm), and the median cold ULIRG 
has stronger absorption than the average starburst (1.36 
versus 0.22). In Figure [3l we plot the apparent silicate 
optical depth as a function of far-infrared spectral slope. 
While the ULIRGs with rg.y > 1 have steeper far-infrared 
spectra than the sour ces with rg y < 1, we fi nd no real 
correl ation (see also llmanishi et al.l (|2007D ; iHao et all 
(|2007f) ). In fact, many sources with very high silicate op- 
tical depths (T9.7 > 2) would also be classified as warm 
ULIRGs. All of the ULIRGs with Scyfert-1 optical spec- 
tra, and many, but not all, of the ULIRGs with Seyfert-2 
optical spectra have rg 7 < 1 and S25/Seo > 0.2. 

The median HII, LINER, and Seyfert ULIRG IRS 
spectra are shown in the left-hand panel of Figurc[4l The 
median and mean values of the 6.2 and 11.3 /tm PAH 
EQW for each of these spectral classes is given in Table 
[TJ As expected, the Hll-like ULIRGs have the largest 
median 6.2 /im PAH EQW (0.28 /im) while the ULIRGs 
classified as Scyferts have the lowest (0.04 /im). The me- 
dian PAH EQW for infrared cold ULIRGs is comparable 
to the median for those classified as Hll-like, while the 
median PAH EQW for infrared warm sources is similar 
to the median for Seyfert-like ULIRGs. There are a num- 
ber of ULIRGs with HH-hke or LINER-hke spectra that 
have unusually low 6.2 /tm PAH EQW (;^0.15 /tm, or 
about 1/4 that found for pure starburst galaxies). Of 
the 22 LINER-like (15 HH-like) ULIRGs in the sample, 
12 (six) have 6.2 /im PAH EQW < 0.15 /tm. Similarly, 
of the 66 cold ULIRGs, 23 have a 6.2 /im PAH EQW 
< 0.15 /tm. 

The right-hand panel of Figure [Jj shows the median HII 
and LINER ULIRG spec tra overplotted with the average 
starburst spectrum from lBrandl et all (|2006f) and the av- 
erage infrared-bright LINER spectrum from lSturm et all 
mm . Although sources with HH-like optical spectra 
have the largest PAH EQWs among ULIRGs, the me- 
dian values are only about 50% of those measured for 
pure starburst galaxies, the same result as found for the 
cold ULIRGs. Hll-like ULIRGs also have stronger sil- 
icate absorption at 9.7 and 18 /im than seen in lower- 
luminosity starbursts (1.52 versus 0.22 for rg.7). Com- 
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pared to infrared-bright LINERs, LINER-like ULIRGs 
have ^ 25% smaller PAH EQWs, and stronger silicate 
absorption (T9.7 = 2.11 for the LINER-like ULIRGs while 
it is only 0.84 for the IR-bright LINERs). 

3.2. PAH equivalent width versus continuum luminosity 

In Figured we plot the 6.2 and 11.3 PAH EQW 
versus 24 /xm rest-frame luminosity. This continuum 
wavelength was chosen because it samples the peak 
emission of warm ('^120 K) thermally-emitting dust, it 
can be measured much more accurately than LIR in 
many cases, it is free from strong emission or absorp- 
tion features, and it allows direct comparison to a va- 
riety of low-redshift samples observed with the MIPS 
24 fim filter and to z 2 samples observed with the 
MIPS 70 fim filter. For reference, we estimate that 
logio(i/L^(24/xm)/LIR) = -1.06 ± 0.47 for cold ULIRGs 
and —0.79 ± 0.38 for warm ULIRGs. These es ti mates 
were made by applying the ISanders fc Mirabell (jl996l ) 
prescription for converting IRAS photometry into LIR 
to the nearest {z < 0.2) ULIRGs in our sample. Figure [5] 
shows a trend of decreasing 6.2 and 11.3 /im PAH EQW 
with increasing rest-frame 24 fim luminosity. There is 
significant scatter in this relation, and ULIRGs with low 
PAH EQW span a wide range in luminosity. Neverthe- 
less, the presence of a correlation is verified by a highly 
significant Spearman correlation coefficient, whether or 
not upper limits arc included in the same way as detec- 
tions in the calculation. 

We have divided the ULIRG sample into six luminosity 
bins chosen to span the full range in 24 fim luminosity 
((0.03 — 8) X 10^^ Lq), and to contain approximately 
equal numbers (17-18) of ULIRGs. The median 6.2 and 
11.3 /im PAH EQWs drop by an order of magnitude be- 
tween the first and last luminosity bins (see Table [2|) . 
The median ULIRG spectrum in each 24 /im luminosity 
class is shown in Figure [6l The median spectra were con- 
structed using all sources in a given luminosity bin, even 
those for which only an upper limit on the PAH EQW 
could be measured. Moving from lowest to highest 24 fim 
luminosity (top to bottom in Figure [B]), the PAH features 
diminish and the spectrum flattens. While the median 
spectrum in the highest luminosity bin has a noticeably 
smaller silicate absorption than the lower-luminosity me- 
dian spectra, there is a large range in the apparent sil- 
icate optical depth among the mo st luminous sources . 
This is consistent with the results of ISpoon et al.l (|2007f ). 
who find a dichotomy in silicate optical depths among 
sources with small 6.2 fim PAH EQW (see also ii3.3p . 

3.3. Silicate absorption versus PAH equivalent width 

Figure [7] shows the silicate optical depth versus 6.2 
fj,m PAH EQW for the U LIRGs in our samp le. Such a 
plot was first presented bv lSpoon et al.l ()2007t ). who dis- 
cussed the presence of two discrete tracks, or branches. 
Figure [7] shows that the horizontal branch is populated 
almost exclusively by infrared warm sources, while the 
diagonal branch is populated by about half of the warm 
sources and virtually all of the cold sources. Examining 
the ULIRGs for which we have optical spectroscopic clas- 
sifications, the diagonal branch harbors about half of the 
Seyfert-like ULIRGs, and nearly all of the Hll-likc and 
LINER-like ULIRGs. The horizontal branch is populated 



mainly by ULIRGs showing Seyfert-like spectra, includ- 
ing all of the Type I sources. ULIRGs with Seyfert Type 
2 optical spectra are found on both branches, spanning 
(nearly) the full range in silicate optical depth. Most of 
the ULIRGs without optical spectroscopic classifications 
lie on the diagonal branch. As noted in §3.11 a significant 
fraction of the HII and LINER-like ULIRGs have very 
low 6.2 fim PAH EQWs (<0.15 fim). Of these, almost 
aU of the LINERs (11/12) and half (3/6) of the HH-like 
ULIRGs have rg.7 > 2.5, implying Ay > 45 mag. 

4. DISCUSSION 

The median 6.2 fim PAH EQWs of cold (or HH-like) 
ULIRGs is about 50% of that seen in lower-luminosity 
starbursts. Warm ULIRGs have 6.2 fim PAH EQWs 
that are up to two orders of magnitude smaller than 
those found in starbursts. There are several possible ex- 
planations for the extremely large range in PAH EQW 
and the overall low PAH EQW in ULIRGs compared to 
starbursts. Some ULIRGs may host a central AGN, and 
the associated soft X-ray an d ultraviolet radiation may 
destr oy PAH molecules (e.g. lAitken fc Rochelll985t IVoitj 
Il992l ). or thermalize the emission. The AGN may also 
heat dust to significant temperatures (>400 K), resulting 
in excess continuum emission which would lower the mea- 
sured 6.2 fim PAH EQW, in efffect diluting the star burst- 
driven PAH emission with respect to the AGN-powered, 
hot dust continuum. For either of these explanations, 
55 ± 13% of the LINER-hke ULIRGs and 40 ± 15% of 
the HH-like ULIRGs in our sample wo uld harbor buried 
AGN. In comparison. iLutz et al . I (| 1999( 1 find the percent- 
age of LINER-hke ULIRGs classified by ISO as AGN to 
be 17± 15% and the percentage of Hll-likc ULIRGS clas- 
sified by ISO as AGN to be 20± 12%. We are thus finding 
a much higher fraction of LINERs and HH-like ULIRGs 
with weak PAH emission, although, given the small sam- 
ple sizes and the fact that the percentage is comparable 
to the statistical errors in the ISO sample, our resu lts 
are formally consistent with those of lLutz et al.l ()1999f ) at 
the 2(7 level. Our results are more comparable to those of 
Iknanishi et al. (2007), who find that 30-50% of HH- and 
LINER-hke ULIRGs contain dominant buried AGN. Of 
the 12 LINER-like and 6 HH-like ULIRGs with low PAH 
EQWs, 9 were observed with the high-resolution mod- 
ules of the IRS (jFarrah et al.ll2007f l. Of these, none have 
detectable [NeV]14.3/im, [NeV]24.3/im, or [OIV]24.9//m 
line emission. Upper limits on the [NeV]14.3/im to 
[NeII]12.8/im line flux ratio range from 0.02 to 0.3, corre- 
sponding to starburst-dom inated ULIRGs wit h less than 
5-25% AGN contribution (jArmus et al.ll2007D . However, 
it is important to note that many of these sources are 
deeply buried, as evidenced by their large silicate optical 
depths (see tj3.3p . As a result, the [NeV]14.3/im emission 
may be highly extinguished. In addition, it is plausible 
that a completely buried starb urst exists in these sources, 
as in NGC 1377 (iRoussel et a l. 2006). 

The 6.2 fim PAH EQW in ULIRGs could also be (in- 
trinsically) smaller than that found in pure starburst 
galaxies if star formation in ULIRGs is accompanied by 
larger quantities of hot dust emission than is typically 
found in lower-luminosity starbursts. The dust would 
then have to "see" more of the stellar radiation and com- 
pete more effectively for UV photons, thus raising its 
average temperature, something that is not unreason- 
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able given the compact nature of the infrared emission 
in ULIRGs. However, if this were the case, then both the 
6.2 and the 11.3 (im PAH EQW would be similarly de- 
pressed. However, among cold ULIRGs, the median 11.3 
fim PAH EQW is 80% of that found in starbursts, while 
the median 6.2 ^m PAH EQW is only 50%. A com- 
bination of dilution and extinction offers a more likely 
explanation. In this scenario, the 6.2 /im PAH EQW 
is diluted by hot dust continuum produced by a central 
ionizing source interior to the star-forming regions. The 
dilution of the 11.3 fj,m PAH EQW is less significant be- 
cause of its proximity to the broad silicate absorption 
feature at 9.7 fiui. The hot dust continuum responsible 
for diluting the 6.2 /xm PAH EQW is heavily extincted 
at 11.3 /im, and therefore does not readily dilute the 11.3 
/im PAH EQW. This is seen in Figured The left-hand 
panels show that the 6.2 /im PAH EQWs of ULIRGs on 
the upper branch deviate strongly from the starbursts. 
In contrast, the right-hand panels indicate that ULIRGs 
have nearly starburst-like 11.3 /im PAH EQWs because 
the strong extinction at this wavelength overwhelms all 
but the strongest contribution from hot dust. 

Finally, extinction is not likely to play a large role 
in the reduced 6.2 /xm PAH EQWs found in ULIRGs 
compared to starbursts. Although ULIRGs on average 
display stronger silicate absorption features than lower- 
luminosity starbursts, larger column densities of dust 
would lead to depressed PAH EQWs only if the PAH 
emission were more extinguished than the continuum. 
Such a situation is unlikely, since the continuum flux is 
emitted by hot dust that must lie close to the ionization 
source. 

In a further attempt to disentangle reduced PAH emis- 
sion from increased hot dust emission as the driving force 
behind the low 6.2 /im PAH EQWs in ULIRGs, we plot 
the 6.2 /im PAH luminosity against the rest- frame 5.5 /im 
luminosity in Figur e [51 Low-luminosity starbursts from 
iBrandl et al.l ()2006f) follow a correlation between these 
two quantities, characterized by logio[L(6.2 /im PAH)] 
= -0.49 + 0.96 X logio[i/L^(5.5/im)]. ULIRGs do not lie 
on the extrapolation of this relation to higher luminosi- 
ties. Instead, ULIRGs appear offset to higher 5.5 /im 
luminosities for their measured 6.2 /im PAH luminosi- 
ties. Not surprisingly, sources with flatter far-infrared 
spectra and those with Seyfert-like optical spectra are, 
as a group, the most displaced towards higher 5.5 /im lu- 
minosity. However, the range of 6.2 /im PAH luminosity 
within the different classes (both optical and infrared) 
is the same. It has been suggested that an evolution- 
ary link exists between ULIRGs and QSOs, such that 
cold ULIRGs evolve into warm ULIRGs o n their way to 
becoming QSOs (c.g. lSanders et allll988a^ . Figure[H]im- 
plies that the transition between cold and warm ULIRGs 
involves primarily an increase in hot dust, rather than 
the suppression of PAH emission. iLutz et all (|l998f l also 
suggested that dilution of the EQW by hot dust was the 
dominant effect among ULIRGs. 

Also plotted in the lower left-hand panel of Figure [8] 
are ~80 PG QSOs at z < 0.5 (Ogle & Antonucci 2007, in 
preparation; Shi et al. 2007, submitted). Most have only 
upper limits for the 6.2 micron PAH EQW, but 15 have 
detected values. The 5.5 /im continuum luminosities of 
ULIRGs and PG QSOs span the same range. The 15 
PG QSOs with detected 6.2 /im PAH features show no 



evidence for a corre lation in Figure [8l This is in contrast 
with the results of ISchweitzer et alT(|2006[ ). who find a 
correlation between the 7.7 /xm PAH luminosity and the 
6 /im continuum luminosity among a sample of 25 PG 
QSOs, including 15 with upper limits on the PAH lumi- 
nosity and 11 with detections. Based on the assumptions 
that the 6 /xm luminosity is due to hot dust heated by 
an AGN, while the PAH luminosity is emitted by star- 
forming regions, they argue that their observed correla- 
tion is evidence for a starburst-AGN connection among 
PG QSOs. However, we find no such correlation when 
considering the 6.2 /im PAH feature. 

While the flatter spectra and lower 6.2 /xm PAH EQWs 
of many of the ULIRGs can readily be explained by an 
excess of hot dust emission, the silicate optical depth, 
at least in some cases, seems inconsistent with a simple 
model for the emission. The 9.7 /xm silicate optical depth 
in ULIRGs varies from —0.2 < Tg.7 < 4, and is on average 
deeper th an in starburst galaxie s. The radiative transfer 
models of lLevenson et alT(|2007[ ) suggest that the geome- 
try of the obscuring material in ULIRGs determines both 
the depth of the silicate feature and the infrared spec- 
tral slope. According to these models, the large appar- 
ent silicate optical depth (rg.y > 1-2) observed in some 
ULIRGs requires an optically and geometrically thick, 
smooth dust component, since a clumpy dust distribu- 
tion can only result in a shallow (rg.y < 1) absorption 
feature. All of the Seyfert-like ULIRGs would then be 
viewed through this clumpy dust distribution, exhibiting 
flat infrared spectra, weak PAH EQWs, and shallow (or 
nonexistent) silicate absorption features. If this simple 
model is correct, a correlation between silicate optical 
depth and spectral slope might be expected. However, 
we find no such correlation (see Figure [3|). In fact, a 
significant number of the sources with the most extreme 
absorption are warm (see Figure [7]). These ULIRGs also 
tend to have low 6.2 /xm PAH EQWs (< 0.15/xm). A 
majority of the sources with large silicate optical depth 
are classified optically as LINERs. 

Among local ULIRGs, the sources with the highest 
rest- frame 24 /xm luminosities have the lowest 6.2 /xm 
PAH EQWs (see Figure [5]). The fit to this correlation 
is shown as the black line in Figure [9l and is given 
by logio(6.2 /xm PAH EQW[/xm]) = (7.71 ± 0.07) -f 
(-0.723 ± 0.006) X logio(i/L^(24/xm)[LQ]). According 
to this fit, ULIRGs with rest-frame 24 /xm luminosi- 
ties above about 10^^ Lq have 6.2 /xm PAH EQWs 
below 0.1-0.2 /xm. However, very luminous infrared 
sources with strong PAH emission (6.2 /xm PAH EQW 
> 0.3-0.5 /xm) have recently been found in high-redshift 
samples that were selected in the sub millimeter, which 
is dominated by cool dust emis s ion (iLutz et all l2005l : 
iMenendez-Delmestre etal\ l2007l : IVahante et al.l boOTt) . 
For example, of the five IRS spectra of submi l limete r 
sources presented by IMenendez-Delmestre et al] (|2007( ) , 
two have the wavelength coverage to measure the 6.2 /xm 
PAH feature. These two submillimeter-selected ULIRGs 
have rest-frame EQWs of approximately 0.6 and 0.8 
/xm, comparable to the strongest PAH emission mea- 
sured for local starburst galaxies. Their positions are 
marked in Figure |9l and they are clearly well above 
the fit to the local ULIRGs. Additional examples have 
been found i n sarn plcs selected at 160 /xm (MIPSJ142824; 
iDesai et al.ll2006| ) or a combination of 24 /im flux and 
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mid- infrared color (|Yan et al.ll2005t ISaiina et all [20071 ). 
These are also plotted in Figure [9l However, sources 
with high rest-frame 24 /zm luminosities and large 6.2 
/Lim PAH EQW are rare in bright (f2 4 > 0.7 mJy) flux- 
limited samples select ed at 24 ^m (e.g. lHouck et al.ir2005t 
IWeedman et al.l[2006[ ): these samples appear dominated 
by AGN. ULIRGs with high rest-frame 24 /xm luminosi- 
ties and strong PAH emission are absent from local sam- 
ples, yet they obviously exist at high rcdshift, and they 
can be found by selecting on color, cold dust emission, or 
by probing to fainter flux levels in the mid-infrared (ob- 
served /24 < 0.3 mJy). At z < 1, the generation of large 
amounts of warm dust seems to require the presence of 
a dominant AGN. The absence of extremely luminous 
starburst-dominatcd ULIRGs at low rcdshift represents 
a real effect, and not simply an observational bias. 

5. SUMMARY AND CONCLUSIONS 

Using a large sample of 107 ULIRGs observed with the 
IRS on board the Spitzer Space Telescope, we investi- 
gate the relationships between the 6.2 and 11.3 /im PAH 
EQWs, continuum luminosity, infrared spectral slope, 
optical spectroscopic classification, and silicate optical 
depth. Our results can be summarized as follows: 

1. There is an extremely large range in detected 6.2 
/j,m PAH EQWs among ULIRGs, spanning more 
than two orders of magnitude from ^^0. 006-0. 9 fj,m. 
ULIRGs classified as starburst-dominated based on 
their far-infrared colors or optical spectra have me- 
dian 6.2 fim PAH EQWs which arc 50% of that 
measured in low-luminosity starburst galaxies. An 
excess of hot dust, which is most prominent in the 
far-infrared warm or Seyfert-like ULIRGs, appears 
to be the cause of the decreased PAH EQW, as 
opposed to extinction or grain destruction. Ap- 
proximately 55% of LINER-like ULIRGs, 40% of 
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small PAH EQWs and far-infrared colors imply 
substantial amounts of hot dust emission. For these 
sources, a model in which a warm source is viewed 
through either a smooth optically thick or clumpy 
dust shell seems inconsistent with the data. 

3. ULIRGs with the largest rest-frame 24 
/im luminosities have the smallest 6.2 fiia 
PAH EQWs. A fit to our sample suggests 
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types are rare at low-redshift, but those that do 
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TABLE 1 

6.2 AND 11.3 fj.M PAH EQWs versus spectral class 







uncorrected 


6.2/xm PAH EQW 


corrected 6.2^tm PAH EQW 


11.3 Atm PAH EQW 


spectral class 


number 


median 


average 


median 


average 


median 


average 






(/im) 


(/im) 




(^m) 


(^tm) 




HII 


15 


0.32 


0.30 ± 0.16 


0.28 


0.27 ± 0.18 


0.39 


0.41 ± 0.20 


LINER 


22 


0.19 


0.27 ± 0.20 


0.11 


0.17 ± 0.14 


0.50 


0.46 ± 0.18 


SI 


6 


0.02 


0.04 ± 0.05 


0.02 


0.04 ± 0.05 


0.03 


0.03 ± 0.03 


S2 


21 


0.05 


0.10 ± 0.11 


0.05 


0.08 ± 0.07 


0.09 


0.16 ± 0.15 


S1+S2 


27 


0.05 


0.09 ± 0.10 


0.04 


0.07 ± 0.07 


0.07 


0.13 ± 0.14 


Cold 


66 


0.30 


0.32 ± 0.18 


0.24 


0.27 ± 0.18 


0.50 


0.52 ± 0.24 


Warm 


41 


0.04 


0.08 ± 0.09 


0.04 


0.06 ± 0.08 


0.07 


0.13 ± 0.13 



Note. — Statistics on the 6.2 and 11.3 fim PAH EQWs of the ULIRGs in our sample, broken down by spectral class (Column 
1). Column 2 gives the number of ULIRGs in each class. Columns 3 and 4 show the median and average 6.2 /.im PAH EQWs, 
respectively, uncorrected for ice absorption. Columns 5 and 6 give the same statistics for the ice-corrected 6.2 fim. PAH EQW (see 
ISpoon et al.l (]2007l ) for details on this correction). Column 7 shows the median and average 11.3 fim PAH EQWs, uncorrected 
for the underlying silicate absorption. All medians and av erages wer e computed by includ i ng sources w ith upper limits a s 
though they were detections. Optical c l assifications are from [Due et_ajjj[199^jjj^ijj£i^ e t al.l (]1997l Il995l) : ICutri et al.l l]1994r ; 
lAllen et aTl (1199111 : lArmus et all (1198911 : [Xleinm an n et al.l (119881) : ISanders et al.l (ll988bD . The far-infrared classifications of cold 
{f I, {25 fim) / f I, {60 fim.) < 0.2) and warm (/^ (25/i.m)//,^ (60/im) > 0.2) are based on the rest-frame 25 and 60 fim values calculated 
as described in 



TABLE 2 

6.2 AND 11.3 fiM PAH EQW versus 24 fiM luminosity 





24 /im 


uncorrected ( 


3.2 fim PAH EQW 


corrected 6.2 ^tm PAH EQW 


11.3 ^irn PAH EQW 


Class 


Luminosity 


median 


average 


median 


average 


median 


average 




(IOI^Lq) 




(/im) 




(/im) 






1 


0.03 - 0.20 


0.43 


0.47 ± 0.16 


0.41 


0.41 ± 0.18 


0.58 


0.61 ± 0.20 


2 


0.20 - 0.30 


0.33 


0.34 ± 0.18 


0.29 


0.30 ± 0.17 


0.56 


0.52 ± 0.24 


3 


0.30 - 0.44 


0.16 


0.23 ± 0.17 


0.10 


0.17 ± 0.13 


0.44 


0.39 ± 0.18 


4 


0.44 - 0.63 


0.14 


0.17 ± 0.11 


0.12 


0.15 ± 0.11 


0.31 


0.30 ± 0.15 


5 


0.63 - 1.15 


0.10 


0.12 ± 0.10 


0.06 


0.07 ± 0.05 


0.19 


0.31 ± 0.40 


6 


1.15 - 8.00 


0.03 


0.04 ± 0.05 


0.02 


0.04 ± 0.05 


0.05 


0.08 ± 0.08 



Note. — Same as Table except the ULIRG sample is broken down into luminosity bins, rather than spectral class. Each 
bin contains 17-18 ULIRGs. 
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Fig. 1. — Left: The median IRS spectra of ULIRGs classified as cold (red) and warm (blue), arbitrari ly normalized for dis play purposes. 
Right: The median spectrum of cold ULIRGs (red) compared to the average starburst spectrum from [Brandl et aTl l|2006l ) (black), both 
normalized at 24 ^m. 
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Fig. 2. — The 6.2 fim (top) and 11.3 fim (bottom) PAH EQW versus spectral slope. The ULIRGs are color-coded by optical spectroscopic 
classification: ULIRGs with Hll-like optical spectra arc plotted as red stars, LINER-likc ULIRGs are shown as green squares, ULIRGs 
with Seyfert 1 optical spectra arc shown as light blue triangles, ULIRGs with Seyfert 2 optical spectra are shown as dark blue upside-down 
triangles, and sources with unknown optical spectroscopic classifications are shown as white circles. The black dots represent the starburst 
galaxies from Brandl et al. (2006). The dotted vertical line represents the division between warm and cold sources adopted for this paper 
(/^ (25)//^ (60) = 0.2), where cold sources lie to the left of the line and warm sources lie to the right. 
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Fig. 3. — The 9.7 fim silicate optical depth versus the mid-infrared spectral slope. The ULIRGs are color-coded by optical spectroscopic 
classification, as in Figure [2] The vertical dotted line is /j^ (25)//^ (60) = 0.2, the division between warm and cold sources adopted for this 
paper. The spectral slope for the ULIRGs was computed f rom t he IRS spectra as described in [[2] and the spectral slope for the starbursts 
was computed from IRAS fluxes taken from lBrandl et aTl II2006I ). 
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Fig. 4. — Left: The colored lines show the median IRS spectra of ULIRGs optically classified as Hll-like (red), LINER-like (green), and 
Seyfert-like (blue). The black line shows the median spectrum of ULIRGs lacking optical classifications. The normali z ation s are arbitrary. 
Right: The median spectrum of Hll-like ULIRGs compared to the average starburst spectrum from IBrandl et al.l II2006I ) (top red and 
black lines respect ively) and the median spectrum of LINER-like ULIRGs compared to the average infrared- bright LI NER spectrum from 
I Sturm et al.l I I2006I ) (bottom green and black lines respectively). For wavelengths longer than 10 /^m, the high-resolution fSturm et al.l I I2006I ) 
spectrum was convolved with a Gaussian profile with a FWHM of 0.2 fim to approximately match the resolution of the low-resolution 
ULIRG spectrum. The comparison spectra were normalized to match the HII- and LINER-like ULIRGs at 24 fim. 
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Fig. 5. — The 6.2 fim (top row) and 11.3 fim (bottom row) PAH EQWs versus 24 rcst-framc luminosity. The ULIRG sample is 
color-coded by warm (/^(25/^m)//^(60/im) > 0.2) and cold (/^(25/^m)//^(60/.tm) < 0.2) mid-infrared spectral slope (left) and by optical 
spectroscopic classification as in Figure |3] (right). The average PAH EQWs of starburst galaxies from IBrandl et al.l 1(2009 ) are indicated 
by solid horizontal linos and the ±1 a dispersions in these averages are shown by dotted horizontal lines. Upper limits (Scr) to the PAH 
EQWs are indicated by downward pointing arrows. The six luminosity classes (1-6; see text and Table [2]l are labeled on the bottom of 
each panel and are alternately shaded for clarity. 
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Fig. 6. — Median IRS ULIRG spectra for the six 24 fim rest-frame luminosity bins shown in Figure [S] arranged such that the luminosity 
increases from top to bottom. The grey bands highlight the positions of the 6.2 and 11.3 fim PAH features. 
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Fig. 7. — The strength of the 9.7 iim siUcate absorption feature versus the 6.2 (left) and 11.3 (right) fim PAH EQWs. The ULIRGs are 
color-coded according to spectral slope and optical classification in the le ft and right panels, r espectively. The colored symbols are identical 
to those in Figure [5] The small black dots represent the starbursts from lBrandl et all I I2006I) . 



14 




13 



V L^(5.5 jim) V L^(5.5 jim) 




V L^(5.5 jim) 

Fig. 8. — Luminosity in the 6.2 fim PAH feature as a function of 5.5 luminosity. In the top panels, the ULIRGs are color-coded 
according to spectral sl ope (left) and optica l classification (right) using the same symbols as in Figure [5] The solid black line is a fit to the 
starbu rsts analyzed by IBrandl et al.l II2006I ) . The bottom panel features expanded x and y axes to show the locations of the IBrandl et al.l 
II2006I ) starbursts (red diamonds) and ~80 z < 0.5 PG QSOs, displayed as cyan squares (detections) and arrows (upper limits). In this 
bottom panel, the ULIRGs are shown as black circles. 
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Fig. 9. — Median (open cireles) and average (filled circles) 6.2 fim PAH EQW versus 24 fim luminosity. Both the medians and the means 
were computed by including sources with PAH EQW upper limits in the same way as detections. The error bars were computed as the 
dispersion in the PAH EQWs in each luminosity bin, and are centered on the average values in this plot. The x-position of each point is the 
median of the 24 fim luminosities of the sources in each luminosity bin. The thick red, blue, and black lines are fits to the cold ULIRGs, 
the warm ULIRGs, and all of the ULIRGs, respectively. Sources with upper limits on the PAH EQW were included in the fits as detections 
with 10% error bars. The horizontal thin and dotted black lines have the same meaning as in Figu re [S] as do the luminosity bin labels 
across the bottom. The green points are the high-rcdshift (z 2) ULIRGs from lSaiina et ahl 1)20071 ) with measured 6.2 ^m PAH EQWs, 
and rest-frame 24 ^m luminosities estimated from the 70 fim MIPS photometry. The large black star is MIPS J142824. 0-1-352619, a bright 
160 ^m sour ce at z = 1.3. The rest-frame 24 fim luminosity for this object was computed by v ia log-log interpolation of the photometry 
presented in lBorvs et al.l II2006I V and the 6.2 fim PAH EQW was taken from lDesai et al.l I I2006I) . 



